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ABSTRACT 

Hydrodynamic simulations of galaxies with active galactic nuclei (AGN) have typically employed 
feedback that is purely local: i.e., an injection of energy to the immediate neighborhood of the black 
hole. We perform GADGET-2 simulations of massive elliptical galaxies with an additional feedback 
component: an observationally calibrated X-ray radiation field which emanates from the black hole and 
heats gas out to large radii from the galaxy center. We find that including the heating and radiation 
pressure associated with this X-ray flux in our simulations enhances the effects which are commonly 
reported from AGN feedback. This new feedback model is twice as effective as traditional feedback 
at suppressing star formation, produces 3 times less star formation in the last 6 Gyr, and modestly 
lowers the final BH mass (30%). It is also significantly more effective than an X-ray background in 
reducing the number of satellite galaxies. 

Subject headings: galaxies: elliptical and lenticular — galaxies: formation — methods: numerical 
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1. INTRODUCTION 

It has been well established that massive galaxies typ- 
ically contain massive black holes (BHs) at their centers, 
and it is widely believed that these black holes, while 
in their so-called "active galactic nuclei" (AGN) phases, 
can have prof ound influence on th e galaxies which 
host them (see iFerrarese fe F ord 2005| for a comprehen- 
sive observational review). This influence is inferred 
from, among other observations, the correlation between 
the BH mass and the b ulge mass of the host galaxy , 
which was first n oticed bvlKormendv fe Richstonen 1995) 
and studied by iMagorrian et al.l (119981) . an d has since 
been repeatedly confirmed (jHaring fc Rixl I2004D . as 
has the somewhat tighter relation between the BH 
mass a nd the host bulge's velocit y dis persion, discov- 
ered b y IFerrarese fc Merrittl (|2000f) and lGebhardt et all 
([2000?) . The inferred BH mass distribution has also been 
convincingly linked to the quasar luminosity function 
(|Yu fc Trema ind [200l. 

Since the work of ISilk fc ReesI (|1998| ). the relationship 
between the AGN and host galaxy has generally been 
characterized as a process of feedback: the galaxy sup- 
plies gas to be accreted by the BH, which emits some 
fraction of the accreted mass as mechanical energy to 
the surroundings, some fraction via the broad-line winds, 
and some fraction as the radio jets. All these processes 
can heat the surrounding gas, and since gas which has 
been heated is less dense and thus accretes more slowly, 
the feedback process is described as "self-regulating" (i.e. 
a negative feedback loop). 



Thus there are two modes by which AGN may 
give energy to their surroundings: the "mechanical" 
mode, where the AGN inflates local gas bubbles which 
expand through t he inter-stellar and inter-galactic 
medi a (ISM/IGM) (jFabian et al.l 120001: IChurazov etaTl 
120021 : iNulsen et all 120051 ): and the "electromagnetic" 
mode, where photons from the AGN accretion re- 
gion for which the neighboring gas has a low optical 
depth escape and directly interact w ith the rest 
of the ISM/IGM (jSazonov et alJ [2001 . Both of 
these modes may then interact with the surround- 
ing galactic gas by either (or rather, b oth) of two 
mech a nisms: energy-b a sed (i .e. heating) (|Silk fc ReesI 
119981: iWvithe fc Loebl I2003D, or momentu m-based 
(i.e. pressure) ~(lFabian. Celotti fc Erlundl [ 2006: 
DcB uhr. Quataert fc Mai l2010f ). ICiotti fc Ostrikerl 
(|2007f ) have included both mechanisms in their ongoing 
work. Thus there are four conceptual components which 
feedback models may take into account: mechanical- 
energy (a.k.a. bubbles), mechanical-momentum (a.k.a. 
winds), electromagnetic-energy (a.k.a. radiative heat- 
ing), and electromagnetic- momentum (a.k.a. radiation 
pressure). The relative importance of these components 
has been much debated among the authors just men- 
tioned. A fifth component, the thin radio jets, emit 
comparable amounts of energy and momentum with 
the other modes, but these intense beams tend to drill 
through the galactic gas, depositing their energy in 
the IGM. This makes jets less relevant as a feedback 
method, unless they precess, or there is a relative 
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velocity between the AGN and its environment, in 
which case significant heating of galac tic gas could occur 
(|Sternberg fc Sokerll200l [Sokenl2009l) . 

There has naturally been much interest in reproducing 
the AGN feedback process via numerical simulations. 
ISpringel et al.1 (|2005h were among the first to do so using 
a high-resolution 3D smoothed-particle hydrodynamics 
(SPH) code. They found that including an accreting 
BH particle which returns energy to the surrounding 
gas reduces the fraction of baryons which form stars, 
and at late times (for massive galaxies) expels most of 
the remaining gas from the host, creating the classic 
"red and dead" elliptical. Since then, many more SPH 
simulations have been done, e xploring various aspect s 
of AGN growth and feedb ack: iPelupessv et al.l (120071) 
iKhalatvan et al.1 (12001) . I Johansson Burkert fe Naabl 
(12009ft. I Johansson. Naab fc Burkertl (|2009f ). and 
iMcCarthv et al.l (|2010l ). to name a few. 

However, even though it has been shown that 
the emitted spectrum of the average AGN power 
EFp, has a strong peak in the X-rays (~ 30 keV: 
ISazonov. Ostriker. fc Sunyaevl 12004) . and moreover the 
ISM/IGM is k nown to be optically thin to h ard X-rays 
in most cases (M orrison fc McCammon1[T98l . all of the 
simulations just mentioned include only the mechanical- 
energy component of feedba ck. Some recen t studi es have 
included other components: Ostriker et al.l (|2010ft exam- 
ine the effect of mechanical-momentum feedback on AGN 
growth in 1 and 2D simulations and find that including 
this component reduces the final BH mass by a factor of 
100; iDeBuhr. Quataert fc Mai (|2010f) insert mechanical- 
momentum feedback into a 3D SPH code (in fact the 
feedback was characterized as a radiation pressure, but 
only applied to the central 0.2 kpc of the galaxies, mak- 
ing it functionally mechanical) to examine its effects on 
the major mergers of disk galaxies, and found that with 
this component BHs self-regulate effectively during the 
mergers, but without dr iving large quantit i es of g as out of 
the galaxy altogether. iCiotti fc Ostriker] (|2007f ). mean- 
while, performed ID simulations with the full electro- 
magnetic mode, considering the heating and radiation 
pressure from the AGN radiation, but without the me- 
chanical mode. They studied the potential cooling flows 
of "recycled" gas from dying massive stars, and found 
that the electromagnetic mode alone was sufficient to 
drive out half of the incoming gas, with the BH only 
accreting 1% of the total (the remainder forming a star- 
burst). However, until now no 3D simulation has been 
perfo rmed with the el e ctrom agnetic mo de inclu ded. 

In lHambrick et~aT1 (|2009[ ) (hereafter [PaperJ) , we ex- 
amined the effect of various ionizing radiation back- 
gro unds on the prop erties of massive elliptical galaxies. 
In lHambrick et all (|2010f ) (hereafter IPaper III ), we ex- 
amined the same backgrounds with respect to the small 
satellite galaxies of those massive ellipticals. Now we 
turn to a different source of ionizing radiation. It is the 
aim of this paper to present a first qualitative look at, 
and hint at quantitative results from, an AGN feedback 
model incorporating an X-ray electromagnetic mode. 

The paper is organized as follows. In §5J we describe 
the numerical methods and parameters of our simula- 
tions, including the various AGN feedback models. In 
$3] we describe the results obtained from those simula- 
tions, in particular the effects of the different feedback 



models on the BHs themselves, on the gas and stars in 
the host galaxy and on the satellite galaxies. Section [4] 
is discussion and conclusion. 

2. SIMULATIONS AND PARAMETERS 



Our simulation code is GADGET-2, as in IPaper III , 
with the only change being the addition of a black hole 
and associated feedback. We use the UV background of 
iFaucher-Giguere et al.l (120091) (wh ich was designated as 
"FGUV" in IPaper II and IPaper Hf ). The simulation in- 
cludes SN thermal feedback but not winds or any stellar 
mass loss. It includes a simple prescription for metal- 
line cooling using cooli ng rates calculated by Cloudy 
(v07.02, last described in lFerland et al.lll998ft . which pre- 
sumes photoionization and collisional equilibrium for the 
gas and metal atoms, and assumes 0.1 solar metallic- 
ity (note that our X-ray feedback code assumes solar 
metallicity, as discussed in §3.4). A self-consistent treat- 
ment of metallicity would increase the amount of cool 
gas available for accretion at later times and thus en- 
hance the differences between our models. Our simula- 
tions do not include optical depth effects or radiative 
transfer, in particular the self-shielding of dense star- 
forming regions from the ionizing backgrounds, although 
those regions would b e optically thin to X-rays regard- 
less. As in Pap er I j . all simulations were performed 
with initially 100 3 each of SPH (i.e. baryon) particles 
and DM particles, with a gravitational softening length 
of 0.25 kpc for the gas and star particles and twice 
that for the dark matter particles. Gas and star par- 
ticles have masses in the range 4 — 7 x 10 5 M©, depend- 
ing on the size of the box; the assumed cosmology is 
(n m Q^,n h /Q, M ,a % ,h) = (0.3,0.7,0.2,0.86,0.65) as in 
IPaper III . 

The BH feedback works as follows. At z — 9 a single 
seed BH particle of mass 1.5 x 10 6 M Q is created at the 
point of minimum potential in the most massive progen- 
itor halo. The seed mass is chosen to roughly follow the 
Magorrian relation for the galaxies at this redshift; the 
AGN behavior is not par ticularly sensitive to the seed 
mass (Hopkins et al. 2006j). The BH grows at a modified 
Bondi-Hoyle-Littleton rate, 
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where Mbh and i>bh are, respectively, the mass and 
speed (with respect to the surrounding gas) of the 
black hole, and p and c s are the density and sound 
speed of the gas in the SPH kernel centered on the BH 
(ISpringel et al.l 120051: IDi Matteo. Springel fc Hernquist 
1 2005 1: iHopkins et al.l [20061 : Uohansson. Naab fc Burkert 
2009). The free dimensionless parameter as, which we 
choose to be 100, represents the fact that the gas den- 
sity at the BH accretion radius is likely to be much 
higher than the average density in the lo cal SPH kernel 
due to the limited resolution of the code (Hopkin s et al.l 
l200l iBooth fc Schavel [2009T ). Future work could be 
done using the more realistic density-d ependent accre- 
tion efficiency of IBooth fc Schavel (|2009[ ), which reduces 
to a = 1 for large accretion radii. The Bondi-Hoyle 
accretion is capped at the Eddington rate, MEdd = 
AirGmpM-Qn/ e r cax ■ When the BH's notional mass as 
calculated from these accretion rates exceeds its true dy- 
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namical mass, the BH particle swallows nearby gas par- 
ticles as needed to make up the difference. We enable 
artificial recentering of the BH particle on the poten- 
tial minimum of the galaxy, since with our modest mass 
resolution dynamical friction is insufficient to keep the 
young BHs centered (though it is at late times, when the 
recentering has no effect). 

We name our three feedback models "BH", "BHX" 
and "BHXRP" . The "BH" model consists only of ther- 
mal energy deposited at each timestep in the SPH ker- 
nel centered on the black hole, in the amount of E = 
€TC r MBHC 2 , where e r = 0.1 is the overal l radiative ef- 
ficiency of the BH (|Yu fc Tremaind 120021) . and e T the 
fraction of the radiation energy output which is assumed 
to be absorbed thermally by the local gas. Our choice of 
er = 0.005, together with our choice of as — 100, is de- 
signed to produce a final (z = 0) BH mass roughly in line 
with the Magorrian relation for these galaxies. Note that 
our er is 10 times smaller than than the value used in 
other GADGET- 2-based simulations (e. g., Springc l et al.l 
120051 : Uohansson. Naab fc B urkert 20Q9|) , which is due to 
slightly different star-formation criteria than the generic 
GADGET-2 (a lower density threshold for star forma- 
tion, a star-formation timescale which is shorter at low 
densities but constant at high densities and never shorter 
than the cooling time, and requiring the gas to fulfill a 
convergent flow criterion and the Bonnor-Ebert criterion 
to form stars), the upshot of which is that our stars form 
at somewhat lower gas densities, effectively lowering the 
BH a ccretion rate and thus its mass. Other simulations 
(e.g.. lOkamoto. Nemmen fc Bowerl 120081 : lOstriker et al.l 
I2010D have used similarly low efficiencies around = 
0.01 — on the other h and, the Overwh elmingly Large 
Simulations (OWLS) (ISchave et al.l[2010h . using a mod- 
ified GADGET-3, find e T = 0.15 to match the Magor- 
rian relation, though also findin g that the value shoul d 
be reduced for lower resolution (Booth fc Schavd 120091 ) . 
Furthermore, observational ev idence suggests a value 
e T « 0.015 (|Moe et al.l I2009D . although since we ne- 
glect important processes like stellar mass loss and chem- 
ical evolution, our value (and others') for er is con- 
strained far more by the simulation than by physics. 
At any rate, we are primarily interested in the rela- 
tive differences between the models with and without 
radiative feedback, which should not be significantly af- 
fected by the exact star-formation prescription used. We 
call the mechanical-energy component "thermal" feed- 
back for short to distinguish it from the X-ray feedback 
below. This energy component is of unspecified origin: if 
it is assumed to be the result of the broad -line wind, one 
shoul d also include the momentum input ( Ostrikcr et al. 
120101) , which has typically been neglected and which will 
not be included here. 

The "BHX" model has the same feedback component 
just discussed, with the same assumed efficiencies, but 
adds another: the electromagnetic-energy component of 
X-ray radiation from the AGN. This radiation is emit- 
ted from the location of the BH particle with a lumi- 
nosity of Lx = ex£ r MBnc 2 , with a bolometric-to-X- 
ray conversion term ex- This luminosity is converted 
to a flux at each gas particle simply by Fx — Lx/^irr 2 , 
with r the distance of the particle from the BH, and 
the flux is converted to a heating rate for the gas us- 



ing Eq. 36-43 in iCiotti fc Ostrikerl (|2007l ) (based on 
iSazonov. Ostriker. fc Sunvaevl 120051) 7 taking only terms 
which are dependent on the ionization parameter £ oc 
Fx ■ These equations are parametrized in terms of bolo- 
metric flux, and implicitly assume ex ~ 0.04 (as well 
as solar metallicity: see §3.4). We do not include ra- 
diative transfer/optical depth effects for this radiation 
(again, the ISM/IGM generally has a low optical depth 
to X-rays). Nor do we include a speed-of- light delay in 
propagation from the AGN across the box, but since our 
box's high-resolution region is only 2 Mpc in radius and 
we are most interested in the central 30 kpc, the effects 
of no delay should be small. 

The "BHXRP" model adds a third component: 
electromagnetic-momentum, the radiation pressure from 
this X-ray flux, by applying to each gas particle a radial 
force away from the BH particle equal to the X-ray heat- 
ing rate at that timestep divided by the speed of light 
c. We neglect the effect of dust, which dominates the 
opacity by several orders of magn itude at temperature s 
where it can exist (T < 10 3 5 K) (jSemenov et al.ll2003lh 
making this component perhaps an underestimate (al- 
though as discussed in §3.3 below, our heating rate may 
be overestimated due to metallicity effects). 

As an additional point of comparison, we run a set of 
simulations designated "BH+X". These have the same 
feedback mechanism as the regular BH runs, but also 
include an X-ray background field (which the others do 
no t). This field is identical to the "FGUV+X" model 
in iPaperlll : see that paper for details. Briefly, the field 
peaks in strength around z — 2 and falls off sharply at 
earlier and later times, to roughly emulate the observed 
quasar background (i.e., the X-ray heat due to galaxies 
other than the one being simulated). This model al- 
lows us to differentiate the effects of X-rays originating 
in local AGN feedback from generic (spatially uniform) 
background radiation. We also compare to a model with 
no AGN activity whatever, "No BH" . 

We use three sets of initial conditions from iPaper ij 
that were designated galaxies/halos A, E, and M, and are 
so designated here as well. These are elliptical galaxies of 
virial mass 1 — 2 x 10 12 M Q . As will be discussed below, 
all the gala x ies ga v e the same qualitative results . As in 
iNaab et all (I2007T) . Uohansson, Naab fc Ostrikerl ()2009f ) 
and IPaper J . we choose galaxy A as our fiducial case, 
noting any differences from the others where relevant. 

Throughout the paper, all distances are ph ysical ex cept 
where no ted, with assumed h — 0.65 as in |Paper"H and 
IPaper ij . 

3. RESULTS 
3.1. Evolution of the Black Hole Mass 

We first examine the effects of these various feedback 
mechanisms on the black hole itself. It is well-known 
that black holes are characterized only by mass, ele ctric 
charge, and angular momentum jBckcnstcin 119981 ). of 
which only the first is relevant here. 

The top panel of Figure [1] shows the mass of the black 
hole for our four feedback models with galaxy A. Also 
shown is a model where the thermal feedback efficiency 
parameter er has been increased by a factor of 2 to 0.01 
(notated as "BH, 2xFB" in the figure). We see that 
although the additional feedback present in the BHX 
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and BHXRP models reduces the final black hole mass by 
15% and 30% respectively, compared to the BH model, 
these changes are significantly less than the factor of 2 
reduction we see in the model with increased thermal 
feedback — and another model with thermal feedback in- 
creased 10 times has a lower BH mass by a factor of 25. 
This is because the thermal feedback directly affects only 
the immediately neighboring gas particles, which are the 
same particles used to calculate the black hole's local 
density and thus its accretion rate, whereas the X-ray 
feedback affects all gas particles, thus making (as we will 
see) its effect on the host galaxy relatively stronger than 
its self-regulation effect, compared to thermal feedback. 

On the other hand, the addition of an X-ray back- 
ground in the BH+X model increases the final BH mass 
slightly (9%), due primarily to more gas being available 
in the merger event at z = 0.7. All the models except 
the enhanced thermal feedback model are consistent with 
the Magorrian relation for this galaxy, which predicts a 
black hole mas s of ~2 x 10 8 M W , d epending on model, 
using the fit of IHaring fc Rixl (|200l) — since the galaxies 
are ellipticals, we approximate "bulge mass" as the stel- 
lar mass within three stellar half-mass radii (~ 18 kpc) of 
the galaxy center. We note here that if the electromag- 
netic feedback (heating and radiation pressure) is used in 
the absence of any mechanical feedback, the black hole 
becomes too large by three orders of magnitude: it seems 
that for our simulations the electromagnetic mode alone 
is insuffi cient for self-regul a tion. This differs from the 
result of iCiotti fc Ostrikerl (|2007| ) , who found effective 
self-regulation using only electromagnetic feedback; we 
attribute the difference to their inclusion of dust opacity, 
which is certainly v ery important in AGN accretion re- 
gions (Scoville 2003), as well as their neglect of infalling 
satellite galaxies and IGM. 

The lower panel of Figure [TJ shows galaxy E, which is 
much the same as galaxy A, except for the major (1:1 
mass ratio) merger event at z 1.5 (10 Gyr ago), which 
drives a huge amount of gas to the BH and causes the 
mass to jump by a factor of ~ 5. Again in this case 
the BH+X model has the largest final mass (by ~1Q%), 
thanks to the X-ray heating leaving more gas available 
to be accreted in the merger event. We also see that 
BHXRP accretes only half the gas of the other models at 
the moment of merger, making up the difference slowly 
over the next several Gyr; it seems that the radiation 
pressure is effective in blowing out the merging gas when 
X-r ay heating and thermal feedback are not, consistent 
with IPeBuhr. Quataert fc Mai (|2010l ). Otherwise, galax- 
ies E and M are much the same as galaxy A, with BHX 
and BHXRP having lower black hole masses by roughly 
15% and 30%. 

Another factor which we do not model here is the cre- 
ation and mergers of black holes in subhalos. For simplic- 
ity, we have created only one BH, in the largest halo at 
z = 9 (which remains the largest halo until the present). 
While we do not expect the lack to change our results 
qualitatively, including AGN in satellites would reduce 
gas mass and star formation in those systems, leaving 
less gas and stellar mass to be accreted by the central 
galaxy in minor and major mergers. The effect would 
be especially strong on systems like galaxies E, with its 
nearly 1:1 merger. 

Figure [2] shows absolute accretion rates for the same 



0.1 0.3 0.6 1 1.52 34 8 



10 8 



10 7 



10 8 



10 7 



~r 



i i n r 




BH 

BHX 

BH+X 

BHXRP 

BH, 2xFB 
J , i i 




5 10 
Lookback Time (Gyr) 



15 



Fig. 1. — The mass of the central black hole for our various 
feedback models, in galaxy A (top) and galaxy E (bottom). Adding 
the X-ray feedback and radiation pressure decreases the final BH 
mass by 30%, while increasing the thermal feedback by a factor of 
2 reduces the final mass by a factor of 2. 

models with galaxy A. The rates are smoothed over 160 
Myr for clarity; the BH accretion shows high variabil- 
ity over timescales as short as a few timesteps (~ 100 
kyr). For the models with lower thermal feedback, the 
rate peaks near lM /yr around 4 > z > 3, and slowly 
declines thereafter to a final value of ~0.003M Q /yr, ex- 
cept for the peak around z « 0.7, which corresponds to a 
moderate (6.5 : 1) merger event for Galaxy A. Of interest 
is the suppression of the accretion spike at z = 3 for the 
BHX model, and all three accretion spikes (at z — 6,4, 3) 
by BHXRP: we see that these feedback modes are effec- 
tive at self-regulation for very high accretion rates. Also 
of note is that for 4 Gyr after the merger BHXRP has 
more accretion than the other models by a factor of 2 
(which will be significant in §3.3). 

In the following subsections we disregard the model 
with enhanced thermal feedback; in all cases its effects 
are the same or weaker than the plain BH model. 

3.2. Impact on Gas 

We now turn to the effects of the various BH models 
on the host galaxy. We expect our X-ray feedback to 
be quite effective at heating gas, and indeed it is. Fig- 
ure [3] shows the temperature distribution of all gas in 
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Fig. 2. — The accretion rate for galaxy A in Mq/yv with the dif- 
ferent feedback models. Accretion peaks at z = 3 and again with 
the major merger at z = 0.7, after which BHXRP has a sustained 
period of higher accretion. The X-ray luminosity is proportional to 
the accretion rate (Lx = ^x^tMbhc 2 ), with lM0/yr correspond- 
ing to 2 X 10 44 erg/s. The sharp spikes at 7 > z > 3 correspond to 
bursts of Eddington-limited accretion. 

the halo A box at z = 3.2 (during the epoch of peak 
star formation) and the present. At z — 3.2, the No BH 
model has the coldest gas, with a mean temperature of 
6.4 x 10 4 K. Adding BH thermal feedback heats the gas 
11% to 7.1 x 10 4 K, and adding an X-ray background 
heats it 24% more to 8.8 x 10 4 K. However the two BHX 
models are more effective by another 9%, giving a mean 
gas temperature of roughly 9.6 x 10 4 K. At z = 0, the 
effect of feedback X-rays is more pronounced compared 
to a background: BHX and BHXRP, at a mean temper- 
ature of 3 x 10 5 K, are 25% hotter than BH+X, which 
in turn is only 14% hotter than BH and No BH. 

Figure 2] is the same as Figure [31 except that it shows 
only virialized gas: gas with a density more than 200 
times the mean baryon density. Here the effect of the 
AGN feedback X-rays is more pronounced, as we would 
expect, since most of the virialized gas is near the central 
galaxy, where the AGN X-rays are strongest. At z = 3.2, 
the three models without X-ray feedback are all within 
10% of each other in mean temperature at roughly 9 x 

10 4 K, while BHX and BHXRP have mean temperatures 
of 1.6 x 10 5 K and 1.3 x 10 5 K respectively, 40-80% higher 
than the other models. At z — the picture is even 
more extreme. BHX and BHXRP have 2.4 x 10 10 M Q 
and 3.0 x 10 10 Mq, respectively, of virial X-ray gas > 

10 5 5 K, more than 4 times the amount that BH+X has, 
and 6 times the amount of the other models. Unlike 
the Warm-Hot Intergalactic Medium (WHIM) , which is 
usually defined as p < 100/5& ([Smith et al.ll2010l ). this gas 
is dense enough to emit significant soft X-rays. 

We estimate the (Bremsstrahlung) X-ray luminosities 
of the various models via 
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Fig. 3. — Temperature distribution of all gas in the box for halo 
A at z = 3.2 (top panel) and 2 = (bottom panel). BHXRP, 
not shown, is nearly identical to BHX; No BH and BH are nearly 
identical in the top panel. The BHX and BHXRP models have a 
10% higher mean temperature than the BH+X model at z = 3.2, 
and 25% higher at z = 0. 

(Evrard 1990) for gas above 2 x 10 6 K over the virial 
volume (a 500 kpc radius) , and find that all models have 
logLx ~ 38 — 39.5, consistent with t he ga s luminosities 
found bv lBoroson. Kim fc Fabbianol ([201 ID for 30 early- 
type galaxies of similar size. The differences between 
BH models is modest: averaged over the three ICs, the 
BH+X, BHX and BHXRP models have increased lumi- 
nosity by factors of 1.6, 2, and 3, respectively, compared 
with the models without X-rays. Adding BH feedback 
also increases the X-ray effective radius: only 33% of the 
total X-ray luminosity for No BH comes from outside the 
central 10 kpc, while 63% does for BH and BH+X, 68% 
for BHXRP and 73% for BHX (again averaging the 3 
ICs). 

3.3. Impact on Stellar Mass 

We naturally expect the hotter gas produced by the 
AGN X-ray feedback to reduce the production of stars. 
The upper panel of Figure [5] shows the star-formation 
rate (SFR) over time for galaxy A, out to a radius 
of 30 kpc. We see only modest differences in the ini- 
tial star-formation peak (see below), but BH+X and 
BHXRP are both effective at suppressing late star forma- 
tion: BHXRP has a lower SFR than the other models by 
0.5 dex for the the last 6 Gyr (z < 0.6), while BH+X sup- 
presses the SFR by 0.8 dex for the last 1.5 Gyr (z < 0.1). 
BHXRP forms 3.4 x 10 9 M of stars in the host galaxy 
after z = 0.5, and BH+X 3.6 x 10 9 M©, which is less 
than half of the roughly 7.2 x 10 9 M formed by the 
other models. The BH+X result is interesting in the con- 



2.5 F 
2.0 

1.5 
1.0 

0.5 

0.0 
2.0 

1.5 
1.0 
0.5 




z = 3.2 



No BH 
BH 

BH+X 

BHX 

BHXRP 



z = 




0.0 

2 3 4 5 6 7 

log T (K) 
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higher mean temperature than the other models by a factor of 1.5, 
and at z = they are hotter by a factor of 2. 
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Fig. 5. — Log of the star-formation rate over time for galaxy A 
with the various feedback models, considering stars within 30 kpc 
(upper panel) and 3 kpc (lower panel) of the galaxy center at z = 0. 
Where the curves drop off the bottom of the plot they go all the 
way to zero. BHXRP and BH+X are both effective at suppressing 
late star formation for the whole galaxy, while in the center adding 
X-ray feedback actually keeps star formation going longer. 



TABLE 1 

Stellar mass results 



text of Pa per 1 where the X-ray background produced a 
strong burst of late star-formation, as gas that had been 
kept hot through the X-ray background peak at z ps 2 
finally cools and flows to the center of the galaxy to form 
stars (the so-called "cooling flow" ) . Here, a cooling flow 
seems to be forming around redshift of 0.3 (see the lower 
panel of Figure [5]) , but the AGN effectively shuts it off. 
The effect with BHXRP, meanwhile, is clearly related to 
its enhanced accretion rate after the merger event: we 
see from Figure [T] that BHXRP accretes less gas during 
a major merger, leaving a residual which forms a few 
central stars (see below), but more importantly powers 
extra feedback for the next several Gyr, suppressing star 
formation at larger radii. 

The lower panel of Figure [5] shows the SFR for galaxy 
A out to a radius of only 3 kpc (the 3D stellar half-mass 
radius for the galaxy is the range 5 — 8 kpc). The results 
here are somewhat different: while all the models with 
BHs have no star formation at the present, BHX and 
BHXRP have nonzero star formation until roughly 1.5 
Gyr ago, compared to 3 Gyr for the plain BH model (and 
2.5 Gyr for BH+X). In terms of mass, BHX forms 1.7 x 
10 s M in the central 3 kpc after z = 0.5, BHXRP 2.8 x 
10 s M Q (the residual from the merger just discussed), 
which is 2 — 3 times more than the 1 x 10 8 M formed 
by BH and BH+X, but still less than a tenth of the No 
BH model (3 .1 x 10 9 M©). This is an expected result: 
as in iPaper J . when X-ray heating is present, more gas 
is available at moderate and low redshift (z < 2) to be 
drawn to the galaxy's center in major mergers and in 
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BHXRP 


5.51 


10.22 


20.21 


0.553 



Note. — Results are for stars within the specified radius of the 
principal halo at z = 0. Masses are 1O 1O M0. The last column is 
the ratio of stellar to total (implied) baryonic mass for the cen- 
tral galaxy. Adding X-rays significantly decreases the total stellar 
mass and galactic stellar mass, while the BHXRP model actually 
increases the mass in the center of the host galaxy, and e*, in 2 of 
3 cases. 

cooling flows, and the AGN feedback, while eventually 
preventing most of this gas from forming stars, cannot 
do so immediately. 

Table Q] shows the stellar mass of the host galaxy for 
the various models and ICs at three radii: 5 kpc, 30 kpc 
and 2 Mpc (which is essentially the whole box). The to- 
tal star formation in the box (column 5) is suppressed 
slightly (5%) by the presence of the AGN, but the X- 
rays, whether from the AGN or the background, reduce 
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Fig. 6. — Star-formation rate over time for the various feedback 
models, considering stars within the whole box (2 Mpc radius) 
of galaxy A (galaxies E and M are similar). BHX and BHXRP 
effectively suppress the SFR peak at 4 > z > 3, while the X-ray 
background reduces star-formation at 1.5 > z > 1. 

the total stellar mass by 15%, three times as much. How- 
ever, the background X-rays and the feedback X-rays 
cause their suppression differently: as we see in Figure 
which shows the star-formation history of the entire box 
for the halo A runs, the two models with X-ray feedback 
suppress the star- formation peak at 4 > z > 3, while 
BH+X is most effective at a somewhat lower redshift, 
1.5 > z > 1. This clearly reflects the relative timing 
of the X-ray flux between the feedback and background 
sources, since the X-ray background peaks in intensity 
at z 2, while the X-rays from feedback peak when the 
BH accretion rate does, at 4 > z > 3. Thus BHX can 
have a strong effect on the SFR peak — recall that in Fig- 
ures [3] &S] we saw that BHX has hotter gas than BH+X 
at z = 3.2 — while later when the local AGN is less active, 
the X-ray background is relatively more effective. This 
difference in timing also produces a substantial effect on 
how the stellar mass is distributed, as we will see in the 
next subsection. 

The stellar mass in the central galaxy (column 4 of 
Table [TJ is more or less as one would expect: thermal 
feedback and X-rays each reduce the stellar mass by a 
modest amount (~10%). To present the same results 
in a different way, the last column of Table [T] shows 
the baryon-conversion efficiency for the host galaxy, in 
this case identified using th e Amiga Halo Finder (AHF: 
iKnollmann fc Knebd 120 09). The value is defined as 
it was in §42 of iPaper HI . e* = M*/M Ha io(^Af /^b) = 
5M*/MHaio for our cosmology, where M* is the halo's 
mass in star particles within 0.1 virial radii of the cen- 
ter and -Mjjaio its total virial mass. As discussed in 
IPaper III , the absolute efficiencies for our models are all 
significant l y high er than the value derived from SDSS by 
IGuo et alJ(l2010h . -0.2, but we arc interested in the dif- 
ferences between the models. The two models with only 
thermal feedback reduce the efficiency by 0.035 from the 
No BH case, while adding the X-ray feedback (or an X- 
ray background in 2 of 3 cases) reduces it by an additional 
0.035. So our most effective model, BHXRP, reduces the 
baryon-conversion efficiency by 7% from No BH. This is 
only one fifth of the reduction which would be necessary 



to bring our simulations in line with IGuo et al.l (|2010l) . 
but still a significant difference. 

We find more unexpected effects in column 3 of Ta- 
ble [TJ which gives the stellar masses for the central 5 
kpc of the galaxy. For galaxy A, the BHXRP model 
has almost 20% more mass than the other AGN models 
(though still 20% less than the No BH case), and 4% 
more than BHX at the 5 kpc radius. Once again this 
turns out to be associated with the SFR peak (though 
a small part comes from the extra late star formation 
discussed above) : we see in the upper panel of Figure [7J 
that BHXRP has significantly more star formation in the 
peak at z w 4. This is clearly related to what we saw 
in the previous subsection: since BHXRP's feedback is 
more effective at self regulation, it accretes less strongly 
at z sa 4 and thus provides less heating to the surround- 
ing gas, which can thus form more stars. I.e., since the 
gas is pushed out by the feedback instead of merely be- 
ing heated, the ability of the gas to form stars is less 
impaired even as the BH reduces its own growth rate. 
This effect is not robust, however: in the smaller galaxy 
M BHXRP has only barely more central mass than the 
other BH models, and in galaxy E it is BHX which has 
the most. But in this case BHX has slightly less stellar 
mass within 30 kpc than BHXRP, so it is only the central 
concentration which has been enhanced. This is related 
to the major merger event which galaxy E experiences: 
the extra residual gas from the early X-ray heating com- 
pared the no-BHX models means more gas is driven to 
the center to form stars during and after the merger, 
which we can seen in the lower panel of Figure [7J Mean- 
while the radiation pressure in BHXRP forces this gas 
back out, as suggested by Fig ure [TJ so it forms stars at 
larger radii, in agreement with lDeBuhr. Quataert fc Mai 

(Hop. 

Figure [8] gives the radial mass distributions for galaxy 
A in the form of circular speeds. As we would ex- 
pect from column 3 of Table [TJ No BH has the high- 
est peak speed and the steepest inner slope, with BH, 
BH+X and BHX lower by about 25% in peak stel- 
lar circular speed and BHXRP in between. We can 
also parametrize the mass distribution with a "Faber- 
Jackson" statistic, a 4 /M^, where a is the (3D) stellar 
velocity dispersion and M+ the stellar mass within some 
radius of the galaxy center (we choose 30 kpc). This 
statistic should be roughly const ant for various elliptical 
galaxies (jFaber fc Jacksonlll976l) , but we find a variation 
of 40% between models, with the No BH case being the 
highest at 0.065 (km/s) 4 /M Q (corresponding to the most 
central concentration), BH+X and BHX the lowest at 
0.047 (km/s) 4 /M Q , and the other models intermediate. 
Similar results are obtained with the other galaxies. All 
these values are lower by a factor of ~ 3 than the observa- 
tional value of 0.22 ( km/s) 4 /Mm deriv e d from the M* — a 
relation compiled bv lRobertson et ahl (|2006f ). but within 
the intrinsic scatter (and quite similar to those authors' 
simulation results for z = 0). 

3.4. Effects beyond the host 

Inspired by the differences between BH+X and BHX 
in their overall sta r- formati on histories (Figure [5]), we 
revisit the topic of lPaper"T3 and examine the low- mass 
slope of the galaxy mass spectrum, as parametrized by 
the a of iSchechteri (|1976l ). which we will call as to dis- 
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Fig. 7. — Star-formation rate over time for the various feedback 
models, considering stars within 5 kpc of the galaxy center for 
galaxy A (top panel) and galaxy E (bottom panel) at z = 0. In 
galaxy A, BHXRP allows significantly more central star formation 
than the other BH models, while the major merger in galaxy E 
causes more star formation for BHX. 
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Fig. 8. — Circular speed for stars (dashed line), gas (dotted line), 
and total (solid line; also includes dark matter and the BH parti- 
cle), for galaxy A with the various feedback models, z = 0. The 
upturn for the total circular speed curves in the innermost 200 pc 
is due to the mass of the BH particle itself (hence it is not present 
in the No BH model). Adding BH feedback lowers the peak value 
and flattens the inner slope, although adding radiation pressure 
reverses those effects somewhat for this galaxy. All BH models 
remove residual gas from the inner 2 kpc. 

ti nguish it from the Bondi accretion parameter as- As 
in lPaper U as is denned for our purposes by n(M)dM oc 

(M/M*) as e- M / M *dM. The results are given in Tabled 
while the No BH and BH models agree within error 
with the values corresponding to their background from 
iPaper III (and BH+X has a modestly steeper slope likely 
related to the suppression of cooling flows), the BHX 



TABLE 2 

SCHECHTER a VALUES AND MAXIMUM 
"BARREN" HALO MASSES WITH VARIOUS 
BACKGROUNDS 



Name 



log /Merit 



No BH 
BH 

BH+X 

BHX 

BHXRP 



-1.09 ±0.04 
-0.96 ±0.14 
-0.84 ±0.12 
-0.44 ±0.13 
-0.55 ±0.08 



8.67 ±0.07 
8.65 ±0.01 
8.76 ± 0.09 
8.95 ± 0.08 
8.88 ±0.03 



Note. — The Schechter-c* values 
for the star particles, maximum "bar- 
ren" halo masses, and overall baryon- 
conversion efficiency with various ion- 
izing background and AGN feedback 
models. X-ray feedback from the AGN 
greatly suppresses small stellar sys- 
tems. 

and BHXRP models indeed have an extreme effect on 
as, flattening it well beyond the obs erved valu e of ~ — 1 
to —0.44 and —0.55, respectively. In iPaperTll we found 
that a increases as any heating is added: from —1.6 with 
no heating, to —1.3 with a UV background, to —1.0 with 
UV and SN feedback, to -0.75 with UV, feedback and 
an X-ray background. As mentioned above, the timing 
of the X-rays seems to explain the difference between 
BH+X and BHX, since the BHX X-rays peak just before 
the epoch of primary star formation in low- mass galax- 
ies (found in IPaper I] to be 3 > z > 2), while the X-ray 
background peaks toward the end of it, and doesn't sig- 
nificantly suppress star formation until z sa 1.5. 

This result suggests that our X-ray feedback might be 
too strong: in fact, our formula converting X-ray flux to 
heating rate has a term (representing the photoioniza- 
tio n heating and line and recombination cooling, Eq. A35 
in iSazonov. Ostriker. fc Sunvaevl 120051) which is linear 
in Z/Zq — i.e. the metallicity as a fraction of solar — 
where solar metallicity is assumed (Sergey Sazonov, pri- 
vate communication). Thus, since the metallicity of the 
ISM could be O.IZq at early times, our heating rate 
co uld be too high by th at factor. Moreover, the work 
of lHui fc Haimarj (|2003f ) suggests that after reionization 
the equilibrium temperature of the IGM is roughly inde- 
pendent of the intensity of the radiation field, depending 
only on its spectrum, so our r~ 2 attenuation factor may 
be less significant than we would naively think. 

The results for /-Meri t are mo re modest. The quan- 
tity /-M G rit is defined in IPaper U : in short, it represents 
the largest halo mass at which halos have an average 
star:DM mass ratio of less than half the global value 
(0.08). (M cr it is the theoretically-calculated virial mass 
whose escape velocity is equal to the sound speed of its 
gas at the epoch when it should be forming stars; the 
effective correction factor / s» 0.75.) Here again , No BH, 
BH, and BH+X agree well with the values that IPaper II 
gives from their ionizing background models. The X-ray 
background models have somewhat higher values, as we 
would expect from their flatter low-mass slopes. 

4. DISCUSSION & CONCLUSIONS 

The effects of X-ray feedback from AGN are manifold. 
We find that X-ray heating and radiation pressure are 
only moderately effective at self-regulation: they reduce 
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the black hole's mass far less than increasing the thermal 
feedback efficiency does, primarily by suppressing bursts 
of Eddington-limited accretion at early times. The model 
with radiation pressure also accretes significantly less gas 
at the time of a major merger, instead accreting it more 
smoothly over the following several Gyr. The X-ray feed- 
back produces a significant reduction in the host galaxy's 
baryon-conversion efficiency compared to a traditional 
feedback model, but only slightly more than a model with 
traditional feedback and an X-ray background. We note 
however that our baryon-conversion efficiency remains 
well above the observationally-derived value; though this 
problem is hardly unique to the present work it remains 
troublesome. On the other hand, less star formation 
would leave more gas available for AGN accretion, which 
would likely enhance the relative effect of electromagnetic 
feedback. 

The enhanced accretion and associated feedback in the 
radiation-pressure model can also sustain a half-decade 
reduction in the star- formation rate of the host galaxy for 
several Gyr after a major merger event, although the gas 
required for this extra feedback leads to more star forma- 
tion in the central regions, which in turn can lead to en- 
hanced central concentration. The AGN X-ray feedback 
also produces a significant mass of virialized, soft- X-ray- 
emitting gas at the present, which the X-ray background 



does not h ave (whe n ordinary AGN thermal feedback is 
present; in iPaper II we found a significant mass in hot, 
dense gas for a model with X-ray background but no 
AGN feedback), which increases both the X-ray lumi- 
nosity and the X-ray half-light radius. In a serendipitous 
final result, we find that this X-ray feedback is also much 
more effective than an X-ray background in suppressing 
small galaxies and thus flattening the low-mass slope of 
the galaxy mass spectrum, due to its feedback's local 
origin making it effective at heating gas some 2 — 3 Gyr 
earlier than the background. 

Since AGN are known to emit X-rays through their 
host galaxies, we view this study as a vital first step 
toward a more complete model of AGN feedback. More- 
over, since the X-ray luminosity of AGN is relatively well 
constrained by observations (though not without intrin- 
sic scatter), there is little need for a new free parameter 
to join the current as and £t (modulo the effects of 
metallicity and dust). Thus we hope that this "new" 
feedback mode will be employed in future SPH simula- 
tions of AGN, since it is both undeniably present and, as 
we have shown, substantial in effect. 

JPO was supported by NSF grant AST 07-07505 and 
NASA grant NNX08AH31G. TN and PHJ acknowledge 
support by the DFG cluster of excellence 'Origin and 
Structure of the Universe'. 
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